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Abstract. While diverse fracture characteristics have been observed in liquid metal 
embrittlement (LME) depending on the solid-liquid metal pairs, the penetration of 
nanometer-thick liquid metal films along the grain boundary has been identified as 
one of the key mechanisms for embrittlement in many classical LME systems, such 
as Al-Ga, Cu-Bi and Ni-Bi. For example, liquid Ga quickly penetrates deep into 
grain boundaries in Al, leading to intergranular fracture under very small stresses. We 
report on a series of molecular dynamics simulations of liquid Ga in contact with an 
Al bicrystal under a constant strain rate. We identify the grain boundary dislocations 
that are nucleated at the grain boundary groove tip and climb down along the grain 
boundary during Ga penetration and characterize their atomic structures based on 
topological method. 
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1. Introduction 



One of the most interesting environmental failure phenomena in polycrystalline metals 
is the degradation of the mechanical properties of ductile metals in the presence of 
an aggressive, wetting liquid metal and an external stress, which is called liquid metal 
embrittlement (LME). This phenomenon presents a challenge in a wide range of material 
processing scenarios (including welding/brazing/soldering, galvanizing, heat treatment, 
hot working, etc.) and is particularly important in nuclear reactors in which liquid 
metals are used as coolants or as spallation targets. 
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A wide range of solid-liquid metal pairs shows LME. Although there are some 
common features of LME systems, several kinds of fracture characteristics have been 
observed in LME depending on the solid-liquid metal couple [1]. Among those, the 
penetration of nanometer-thick liquid metal films along the grain boundary has been 
identified as one of the key mechanisms for embrittlement in many classical LME systems 
(e.g., Al-Ga, Cu-Bi and Ni-Bi). The Al-Ga couple is probably the most widely studied 
among the systems known for the formation of penetrating intergranular films. It 
exhibits strikingly rapid Ga penetration along grain boundaries, resulting in a very rapid, 
dramatic loss of Al ductility. Many experimental studies show that liquid Ga penetrates 
into Al grain boundaries at a remarkable rate (^several /xm/s at room temperature), 
leading to a distinct channel morphology [21 |3l HI El [6l [7] . The penetration of liquid Ga 
along the grain boundaries produces wetting layers with thicknesses reported to range 
from a few monolayers [2] to several hundred nanometers |H El |6] , even in the absence 
of an applied load. Interestingly, the rate of propagation of such liquid layers is strongly 
influenced by even very small stresses O ^ . 

Several models have been proposed to explain the driving forces and atomic 
mechanisms by which the liquid phase penetrates quickly along grain boundaries [H 
[9l [To]. However, many fundamental questions remain unsettled for the Al-Ga couple. 
For example, none of these approaches successfully explains the effects of stress on liquid 
film penetration. Recently, Nam and Srolovitz [TT] examined LME through the lens of 
molecular dynamics (MD) simulations of an Al bicrystal in contact with liquid Ga and 
observed how Ga penetrates along the grain boundaries with and without an applied 
stress. The simulation study proposed a mechanism of LME and clarified how it is 
affected by applied stresses. The interplay of stress and Ga penetration leads to the 
nucleation of a train of dislocations on the grain boundary below the liquid groove root 
which climbs down the grain boundary at a nearly constant rate. 

In this study, we carry out MD simulations of stress-facilitated Ga penetration 
under a constant strain rate condition. We further investigate on the characteristics 
of the grain boundary dislocations that are consistently observed in our series of 
studies [m [121 US]: We analyze atomic structures of the grain boundary dislocations 
in more detail, based on the topological approach by Hirth and Pond [H]. Our results 
show that the characteristics of climbing grain boundary dislocations would depend 
intrinsically on the grain boundary type/structure. 

2. Simulations Procedure 

All of the simulations were performed on a 3-dimensional Al bicrystal sample in contact 
with liquid Ga, as shown schematically in Fig. [H Liquid Ga was initially presaturated 
with Al (the solubility limit as a function of temperature for these interatomic potentials 
was determined in Ref. [I2]). The atomic interactions in the Al-Ga system were 
described by semiempirical embedded-atom method (EAM) potentials developed in the 
previous study [12] . We imposed periodic boundary conditions in the x- and y-directions. 



3 



liquid Ga 
solid M 
^ (JB ^ 

Fixed layer ^ 




Figure 1. The simulation cell contains two Al grains 320, 000 atoms) in contact 
with liquid Ga 40,000 atoms). Panels (a) and (b) show the atomic structures of 
(a) the S5 36.9 °(301)/[010] symmetric tih boundary and (b) the S17 28.1 °(401)/[010] 
symmetric tilt boundary, respectively. 



fixed several atomic layers at the bottom of the bicrystal (to prevent grain rotation so 
that it mimics the surface region of a macroscopic sample) and left the top surface free 
(i.e., there is a vacuum above the liquid) [12]. (For more details, see Ref. [T2].) 

In order to investigate the stress-facilitated Ga penetration, the MD simulations 
were performed under a constant strain rate condition [NVT ensemble) of e^x = ^o + ^xx, 
where the strain rate ixx was controlled by stretching the simulation cell every 1 ps at 
a constant rate in the x-direction (eo is a initial strain). We performed the simulations 
at a constant strain rate rather than fixed grip or fixed load conditions in order to 
complement the artificial grain size effect. (More discussion is provided in the next 
section.) The simulations were performed at 600 K with a strain rate condition of 
eo = 0.01 and = 1 x lO^s"^ 

A couple of different grain boundary types were examined including 
S5 36.9°(301)/[010] symmetric tilt boundaries, S17 28.1 °(401)/[010] symmetric tilt 
boundaries, and low angle tilt boundaries. Both S5 and E17 tilt grain boundaries have 
relatively high grain boundary energies relative to the low angle tilt grain boundaries. 

The MD simulations were performed using the Large-scale Atomic/Molecular 
Massively Parallel Simulator (LAMMPS) code [IHllIE]. The equations of motion were 
integrated using the velocity Verlet method. Total simulation time was at least 60 ns 
(~ 2 X lO^At, where the time step At =2.5 fs). 

3. Result and Discussion 

We performed a series of MD simulations of an Al bicrystal in contact with liquid Ga 
and examined the initial stage of Ga penetration along the Al grain boundaries at 600 K 
under a constant strain rate condition of eo = 0.01 and ixx = 1 x lO^s^^. The imposed 
strain rate condition corresponds to change of strain from 1.0% to 1.6% for our 60 ns 
simulation time. Figure l^a) shows the Ga concentration and stress profiles at 10 to 
60 ns for simulations of the ZI5 36.9 °(301)/[010] symmetric tilt boundary. In Fig. |2|^b), 
the Ga penetration depth vs. time is shown, where the penetration depth is arbitrarily 
defined as the distance from the surface (in the z-direction) where the Ga content is 
equivalent to 0.3 monolayers (a typical Ga content profile was shown in Ref. [12]). Our 
MD simulations under a constant strain rate condition show similar penetration behavior 
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Figure 2. (a) Ga concentration (left panel) and stress u^x (right panel) profiles at 
iO to 60 ns for simulations at 600 K at a constant strain rates of ^ 1 x lO^s^^. The 
Ga penetration depth vs. time is shown in (b). 



of liquid Ga along the Al grain boundary as under fixed grip conditions [12] . Application 
of a constant strain rate significantly promotes liquid metal penetration by forming a 
set of climbing grain boundary dislocations, resulting in fixed rate penetration mode. 
However, we could observe nucleation and climbing of more grain boundary dislocations 
than in the previous MD simulations under fixed grip conditions [12]. Interestingly, 
grain boundary dislocations are generated continuously all over the simulation time. 

In fact, implement of appropriate stress condition is a central issue in our simulation 
of the stress-facilitated Ga penetration. For example, difference between constant stress 
(NPT) and constant strain (NVT) ensemble simulations can be quite significant in 
this situation. If we implemented a constant stress algorithm (fixed load condition), 
the grain boundary traction within the simulation cell would never drop all over the 
simulation time (with no grain size effect). However, the grain boundary would pull 
apart as liquid Ga penetrated only a few tens of nanometers. This is clearly unphysical, 
since our simulation cell is meant only to mimic a very thin layer near the surface of 
a very thick polycrystalline sample, where most of the load is carried by the solid, far 
from the surface. On the other hand, if we performed simulations under a constant 
strain (fixed grip) condition, the applied stress effect would fade out in relatively short 
simulation time since the initial stress within the fixed system would be quickly relaxed 
by even very small amount of Ga penetration along the grain boundary. This unintended 
grain size effect is unavoidable because the inter-distance of the grain boundaries (under 
the current simulation cell size and periodic boundary condition) is too small compared 
to macroscopic size. In order to reflect real grain size sample, we might carry out full- 
size multiscale simulations, e.g, where the displacement field for a macroscopic sample 
is calculated by using a continuum approach and the atomic process that occurs in 
the embedded region (where the grain boundary intersects the solid-liquid interface) 
is simulated by molecular dynamics. While this could, in principle, be done, the 
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Figure 3. (Color) (a) Comparison of th.e MD simulation (right) and linear elastic 
stress field (left) (b) Atomic scale image of Ga penetration along E5 (301)/[010] 
symmetric tilt boundary in an Al bicrystal. The atoms shown in blue represent Al 
atoms and those in red are Ga. (c) Enlarged view of dislocation core region (Burgers 
circuit). Atomic structure was quenched in order to remove thermal noise. The Grey 
and yellow circles represents rows of atoms with positions in alternating (020) planes. 



computational cost and complexity is too expensive to implement. Instead, Namilae 
et al |T7j applied artificial crack tip opening loads on Al bicrystals (with a pre-existing 
notch) in the presence of hquid Ga, where regions adjacent to the wedge were subject 
to a constant displacement rate. In this way, they studied the effect of Ga on the 
crack opening of Al bicrystals using atomistic simulations. Here, we performed MD 
simulations under constant strain rate conditions as a compromise. 

Clearly, the grain boundary traction in real sample would be relieved as Ga 
penetrates along the grain boundary. However, the rate of stress relieving near grain 
boundaries in an macroscopic polycrystalline sample would be much slower (by orders of 
magnitude off) than we could expect in a constant strain simulation of nanometer . In 
this respect, application of appropriative constant strain rate starting from a finite strain 
value is expected to compensate the size discrepancy in a constant strain simulation. 
In actual situation, the strain rate of embedded region would be a time dependent 
variable coupled with both stress redistribution by Ga penetration at the embedded 
grain boundary region and accompanying relaxation within the surrounding large grains 
far from the grain boundary. However, based on the linear Ga penetration kinetics of 
our previous simulation, we expect that the strain rate near the grain boundary would 
be approximately constant within the initial stage of Ga penetration and treated as a 
magic number parameter for simplicity. With the application of appropriate constant 
strain rates, the nucleation of grain boundary dislocations was dramatically increased 
(up to more than 5 dislocations within our 60 ns simulation time), while the stresses 
within the system remained steady by balancing the displacements introduced by Ga 
penetration into the grain boundary and the increase of the simulation cell size. 

Next, in order to characterize the grain boundary dislocations, we examined the 
atomic configuration near the stress concentration patterns. Figure [3](b) shows the 
atomic scale image of Ga penetration along E5 symmetric tilt boundary captured from 
the box region indicated in Fig. [3]^a). The atoms shown in blue represent Al atoms and 
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those in red are Ga. After quenching the atomic structure in order to remove thermal 
noise, atomic scale image was further enlarged as shown in Fig. |3]^c). Examination of the 
atomic structure of the grain boundary shows the existence of an interfacial dislocation 
at the location of the center of this stress pattern. We could accurately determine the 
Burgers vector from atomic scale images of the dislocation by applying the topological 
theory of linear defects at interfaces suggested by Hirth and Pond [l4j. The Burgers 
circuit was indicated in Fig. Et^b) and the analysis approach is shown in Fig. |3]^c) for 
the S5 symmetric tilt boundary examined in this study. Using this approach, we found 
that this LME dislocation has Burgers vector b given by the following relation: 

b = tx-t, = ^[30lU (1) 

where t\ and are step vectors in the crystal systems of the grain A on the left and 
the grain fi on the right, respectively. 

The Burgers vector is indeed the displacement shift complete (DSC) vector of 
b = ^[310] A in the crystal system of the grain A. The component of Burgers vector 
normal to grain boundary is now bn = \b\ = — 1-28 A, where oq is the cubic lattice 
parameter. We confirm this measurement by analytically calculating the stress field 
associated with such dislocations [18] using the measured bicrystal elastic properties and 
comparing it with the stress field determined in the simulation, as shown in Fig. |3]^a). 
The excellent correspondence confirms that the normal component of the Burgers vector 
bn is ~ 1.28 A. 

A series of [010] symmetric tilt boundaries was examined in the present simulations: 
S5 [010] symmetric tilt boundaries, S17 [010] symmetric tilt boundaries, and low 
angle tilt boundaries. We found that the Ga penetration behavior is sensitive to grain 
boundary type and structure: no grain boundary wetting was observed within our 60 ns 
simulations for the low angle (~ 10 °) grain boundaries, while the S5 [010] and S17 [010] 
tilt boundaries showed remarkable Ga penetration rates with formation and climbing 
of LME dislocations. We also analyzed atomic structure of LME dislocation for the 
S17 [010] tilt boundaries in the same procedure as described above and found that 
the Burgers vector is also the DSC vector with the normal component of the Burgers 
vector bn = \b\ = ~ 0.98 A. The Burgers vector analysis implies that characteristics 
of LME dislocation should be dependent on grain boundary crystallography and this 
explains why Ga penetration kinetics is sensitive to grain boundary type/structure. 

Similar to the previous reports [HI [T21 US], our simulations demonstrate that 
application of a constant strain rate significantly promotes liquid metal penetration 
along grain boundaries, resulting in a change from a diffusive to fixed rate penetration 
mode. First, Ga diffuses down the grain boundary in Al below the liquid groove root 
and causes stresses large enough to nucleate a dislocation in the grain boundary. The 
first dislocation "climbs" down by stress-enhanced Ga hoping across the dislocation 
core, leaving a tail of Ga behind. This Ga hopping leads to a constant dislocation 
climb rate that is applied stress- independent. Once the dislocation moves far enough 
from the groove root, another dislocation is nucleated. It too climbs down the grain 
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boundary at the same rate, resulting in a uniform spacing of climbing dislocations. With 
Ga at the grain boundary, applied strains enhance the grain boundary opening and in 
turn more Ga is inserted from the liquid groove into the grain boundary to relieve the 
residual stress (i.e., Ga layer thickening process). The Ga penetration rate mirrors the 
dislocation climb rate and hence is time independent. Here, the applied strain rate 
provides stress field to aid the nucleation of dislocations at the grain boundary and 
keep the grain boundary open to allow sufficient Ga transport enough to move with the 
dislocation. 

4. Summary 

We performed MD simulations of an Al bicrystal in contact with liquid Ga and examined 
the penetration of Ga along the Al grain boundaries under a constant strain rate 
condition. Consistent with the previous reports [HI [121 [13] , we could observe that LME 
dislocations are nucleated and climb down along grain boundary, resulting in fixed rate 
penetration mode. With the application of a constant strain rate, the nucleation of 
grain boundary dislocations was dramatically increased (up to more than 5 dislocations 
within our 60 ns simulation time), while the stress within the system remained steady 
by balancing the displacements introduced by Ga penetration into the grain boundary 
and the increase of the simulation cell size. We also analyzed the grain boundary 
dislocations based on topological approach and characterized it in the framework of 
interfacial structure. Our analysis on several different grain boundaries show that the 
characteristics of climbing grain boundary dislocations depend on the grain boundary 
type/structure and this implies that Ga penetration rate could be sensitive to grain 
boundary type/structure due to the different topological properties of grain boundary 
dislocation. 
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